Block model for the XY-type Landau-Ginzburg- Wilson Hamiltonian with an 

inhomogeneous temperature 
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The phase fluctuation near the saddle point solution of the XY-type Landau-Ginzburg- Wilson 
Hamiltonian with random temperature is studied. For the modes with lowest eigenvalue, the systems 
is self-organized into blocks, which are coupled as a XY model with random bond. The couplings 
obtained in this way agree with those by domain wall method. 
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I. INTRODUCTION 

In recent years more and more experiments show lo- 
cally ordered regions (LOR). Localized Bose-Einstein 
condensation above the global superfhhd transition tem- 
perature is revealed in superfluid transition of He in 
silica gel p], 0]. It is well-known that for some granu- 
lar superconductors, on the insulating side of the su- 
perconductor insulator transition, each grain is sepa- 
rately and independently superconducting while a trans- 
port measurement shows the film to be insulating [3|. 
In recent experiment on amorphous NbN films, pseudo- 
gapped state due to the locally superconducting islands 
is discovered Q. The nucleation of pairing gaps in 
nanoscale regions above T c is also found in the high- 
T c supercondcutor BiiSr^GaCu^O^s [E]. In addition, 
a local metallic state is observed in globally insulating 
La\.2iSri/j(,Mn20j well above the metal-insulator tran- 
sition [6]. The existence of Ferromagnetic region in the 
paraphase of Lai- x Ba x MnOs is discovered [7J. 

To understand the relation between LOR and the 
phase transition in a general way, the saddle point equa- 
tion of Landau-Ginzburg Hamiltonian with random tem- 
perature is solved recently Q . LOR is explicitly shown 
in these solutions. Moreover it is found that there exist 
many excited solutions, which minimize the Hamiltonian 
locally in the configuration space. These solutions can 
be described by the block model (in the following we 
call it B model), in which the system is self-organized 
into blocks. These blocks behave like superspins and are 
coupled with their neighbors. In reference pxil ]. a gen- 
eral method to calculate the couplings between adjoined 
blocks is proposed. This method is based on the free 
energy increasing of domain wall (DW) between the ad- 
joined blocks. So we call this method DW method. 

However for the systems with continuous order param- 
eter, such as superfluid, superconductor, the continuous 
phase fluctuation about the saddle point solution should 
be taken into account. In DW method this kind of fluc- 
tuation is absent. In this paper we study the XY-type 
landau- Ginzburg- Wilson (LGW) Hamiltonian with ran- 
dom temperature. We propose a Gaussian approxima- 
tion to study the continuous phase fluctuation near the 
saddle point solutions. We will show that 



(1) . The system is still be organized into blocks given 
by the DW method. 

(2) . For the modes with lower eigenvalues, the cou- 
plings between the blocks are XY-type like a Josephson 
junction array and the couplings are approximately equal 
to those given by the DW method. The blocks are cou- 
pled like XY-type spins. 

Our paper is arranged as follows. In section II, 
the model of XY-type landau- Ginzburg- Wilson (LGW) 
Hamiltonian with random temperature is given. In sec- 
tion III, one dimensional case is discussed. In section IV, 
two dimensional case is discussed. Section V is a sum- 
mary. 



II. THE MODEL 



We consider the XY-type landau- Ginzburg- Wilson 
(LGW) Hamiltonian with random temperature 

H = J dv{\\ V 0(r)| 2 + \t{v)tf{v) + (1) 



where 



iv^r+iv^i 2 , 



( 2 ) 

and t(r) — t + t(r), and t,t(r) are the average reduced 
temperature and the random part caused by the disorder 
respectively. The parameters are scaled according to the 
references jl, [l(| ■ 

The saddle point equations are given by 

- V 2 ^(r) + [t(r) + 2 (r)]^(r) - 0, (3) 

-V 2 ^(r) + [i(r)+0 2 (r)]^(r)=O. (4) 

Through this paper we assume the saddle point solu- 
tions are along (p x direction, 



0jM 4> v = 0. 



(5) 



where v is used to label the excited states. 

Substituting the saddle point solution into Eq. (J]), 
one get the free energy |Tl| 



F V =H{{^}) = - J dv\{^)Hv), 



(6) 
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for the i/th solution. 

If we assume the solution is along <p x direction and 
ignore the fluctuation in (f> y direction, the problem is 
reduced to be Ising-type. It has been shown that in 
that case the system is self-organized into blocks and 
the blocks are coupled like Ising-spins |9J . The couplings 
between blocks can be obtained by calculating the free 
energy increase due to the domain wall [ldj . If ((j> x ,<fiy) 
is regarded as a complex parameter, letting <j> y — is a 
constraint that only the phase of and it is allowed. If 
we take the fluctuation of <f> y is into account, the phase 
fluctuation becomes continuous. 

In order to write the Hamiltonian in terms of the am- 
plitude and phase of the order parameter, we introduce 

4> x — $cos0, = <£sin0. (7) 

Then LGW Hamiltonian becomes 

H = J dv{ l -[\ V ^ + ^\ V 9\ 2 ]+ l -t{r)^ + \^}. (8) 

In this form, we can see that the free energy increase 
induced by the variation of phase is proportioned to the 
square of amplitude. 



III. ONE-DIMENSIONAL CASE 
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FIG. 1: (Color on line) (a) The ground state and two excited 
state solutions with temperature field with tw = —30, tb = 
10. (b) The first 15 eigenvalues of Eq. (JTTJ) . (c) 4 typical 
eigen functions, (d) The phase variations corresponding to 
the eigenfunctions in (c). The straight line is &%. 



A. The method of domain wall 

As an example, we first consider a system with size 
being 14 consisting of 7 wells and 7 barriers, i.e. the 
temperature field is given by 

■«>-{£ V?4f s V < 9 » 

where i = 1, 2, 7. Here we let the spatial coordinate be 
£ to distinguish from the directions of order parameter. 

Obviously this system is periodic with period of 2. The 
ground state solution <fi x ^ for the temperature field with 
t w = — 30, ti, = 10 is shown in the Fig. (la). The saddle 
point equation is solved by finite-difference method with 
step h = 0.025 [12]. The saddle point solution is approx- 
imately equal to \J —t w at the centers of wells and decays 
to very small in the barriers [8j. Using the DW method, 
we can show that there are 7 elementary blocks and can 
calculate the couplings between the adjoined blocks. For 
example, letting the initial value of (f> x be negative in the 
first well and positive at other positions, we will get the 
excited solutions (f>x , which is also shown in the Fig. 
(la). The spatial range of the first block is given by 

<f£\t) < 0. Explicitly it is 0.5 < £ < 2.5. 

Similarly, letting the initial value of <p x be negative 
in the second well and positive at other positions, we 

(2) 

will get the excited solutions <fi x . The spatial range of 
the second block is 2.5 < £ < 4.5. We do not show it 
in figure since its shape is the same as (jJiP and can be 



obtained by shifting <p£' in £-axis by 2. Therefore the 
spatial ranges of the 7 elementary blocks are given by 
2(i - 1) + 0.5 < £ < 2i + 0.5, where satisfy (0 < 0, for 
i= 1,2,--- ,9. 

Here we give a general method to get the spatial range 
of the elementary block. For the general disordered cases, 
the domain wall of elementary block can obtained by the 
method of opening windows fib} . The spatial range of 
this block is that surrounded by the domain wall. 

Letting the initial value of <j) x be negative in both the 
first and second well and positive at other positions, we 

(12) 

will get the excited solutions tfi x . Substituting the 
solutions <ffi ,$x\<t>x ,<t>x^ into Eq. ([6]), we get the 
free energies F , Fi,F 2 ,Fi 2 , and the free energy increases 
A = F 1 - F , h = F 2 - F and f 12 = - F . Then 
the couplings between blocks is given by [Toj] 

K[ 2 ] = (fi + h - /ia)/2. (10) 

Here we use "(D)" to denote the DW method. In this 
way we get the couplings between adjoined blocks. Then 
the free energy of z^th state is given by 

fu = - Yl #M+ifo°-i+i - !)/ 2 ( n ) 

i 

where oi is the sign of the ith block in the z^th solution. 

The couplings between adjoined blocks with t w — —30 
at different tb are given in table 1. 



3 



B. Gauss approximation: expansion around the 
ground state 

Considering the symmetry of the system, the effective 
Hamiltonian should be XY type rather than Ising type. 
Simply let c^, <jj be unit two dimensional vectors, i.e. 



o = (<r x , a v ), 



°l + ^ = 1 



(12) 



Eq. (|TTj) becomes a XY model. DW method does not 
contradict to this conclusion and it only provided us the 
excited states with a x = ±1,(7- = 0. 

To study the cases with continuous a, we expand the 
GLW Hamiltonian near the ground state saddle point 
solution. Let 

4> x = 4> x - <fi x 0) , 4> v = 4> y - 4 0) = ^ . ( 13 ) 

and wc consider the Gauss approximation that 

H^F + 5H X + 5H y (14) 
where Fq is the free energy of ground state solution and 

5H X = jf" \[{ d -^f + [t + 3(*£") a ]#]de (15) 



SH, 



14 



1, 



«-/ ^[(^) 2 + [i + (4 0) ) 2 ]^]^ (16) 



where the quartic terms are omitted. 

Because the saddle point solution is assumed along (f> x 
direction, the fluctuation of <fi x is the amplitude fluctu- 
ation and that of <fi y is just the phase fluctuation. The 
eigenmodes of <fi y satisfies the following equation 

-^ + [i + (4 0) ) 2 K = A„^„ (17) 

where A„ is eigenvalues and ip„ are the eigenfunctions. 
Then we have 



where 



14 

"„ = / 4>y(p n d^. 

lo 



(18) 



(19) 



The contribution of phase fluctuation to the partition 
function is given by 



(20) 



We discretize the equation (fT7|) with grid of step 0.025 
and solve it by LAPACK, which is a package to deal 
matrix. The first 15 eigenvalues are shown in Fig. 1(b) 
and some eigenfunctions are shown in Fig. 1(c). 



The modes of <j) y with lower eigenvalues give the main 
contribution beyond the saddle point solution. As shown 
in Fig. (lb), where U = -30,4 = 10, the first 7 
eigenvalues are remarkably lower than other eigenvalues. 
Then they possess much bigger thermodynamic ampli- 
tudes than other modes according to Eq. (f20|) . 

In order to understand the first 7 eigenfunctions more 
clearly, we introduce 



tan9„ = <p nl 



(21) 



and show 0i, 63, 67, 6s in Fig- (Id). This function show 
the phase variation of the eigenfunctions. As one can see 
in Fig. (Id), 0i is a constant. This is because that the 
first eigenfunction ipi has eigenvalue Ai = and satisfies 



14 



M\2 



(22) 



This can be shown by comparing the equation (3) and 
(IT71) . This eigenfunction corresponds to a global rota- 
tion, so its eigenvalue is zero. Ai = corresponds to the 
infrared divergence. 

Observing 63,67, one can see that the variation of 
phase in the wells are obviously smaller than in bar- 
riers. 82,64,65,66 also have the this feature. From 
Eq.©, we can see that the free energy increase related 
to the phase fluctuation is proportional to the square of 
ground state saddle point solution. In the barriers, the 
saddle point solution is much smaller than in the wells. 
Therefore phase fluctuation in the barriers induced small 
energy increase. The eighth eigenfuction does not have 
such a feature and its energy is remarkably higher than 
the first 7 modes. This means that for the first 7 modes, 
each block can be regarded as a unit. The phase varia- 
tion inside the block can be ignored and only the phase 
difference between blocks are concerned. Therefore we 
introduce the block functions 



,(o) 



(0; 



0; 



other 



<0, 
cases. 



(23) 



H) 

where the spatial range given by (fix < is explicitly 
given by 2(i — 1) + 0.5 < £ < 2i + 0.5 as mentioned in 
subsection III A. 

Then we assume that 



(24) 



;.=i 



where only one phase is assigned to each block, then Eq. 
(US) becomes 



where 



7 7 

SH V X! ^1 J :i ' <l J 

i=l j=l 



I 7 

n=l 



(25) 



(26) 
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with 

/■14 

Ai, n = / *i(£Vn(0<£- ( 27 ) 

Jo 

In this effective Hamiltonian, we only take the first 7 
modes into account. 

On one hand the effective Hamiltonian Eq. (|2"5|) can 
be given by the expansion of the following XY model 
approximately 

SH y « - £ 4 G) (cos(^ - 0;) - l)/2 (28) 
for 6»i, 6»j < 1 with 

</m+i = Ji+i,i = -K^lJA, 
Ji,i+i = J i+ 2,i = -K[%/A (29) 

where "(G)" is used to denote the method of Gauss ap- 
proximation. This Hamiltonian is consistent with Eq. 
(fTTj) obtained by DW method. If we regard Eq. (fTT]) is 
a special form of Eq. (f28| with 6i = 0, 7r, it should have 

^ii+i = -f^+i- We investigate the cases with different 
tf, with fixed t w — —30. The numerical results for matrix 
elements of Jij and the couplings K^h by DW method 
are shown in Table 1. As shown in table 1 and Eq. (|29|) 
, i+i ~ Kii+i i s satisfied even for tb — in an error 
less than 20%. For higher tb, two methods agree with 
each other very well. At tb = 40.0, the relative difference 
between and is less than 10~ 5 . In addi- 

tion the couplings between next nearest neighbors are 
much smaller than the that between nearest neighbors, 
i.e. Jj,j+2 Ji,i+i- This indicates the approximation of 
nearest neighbors is good enough. 

Moreover Eq. (|2"51) is the well-known Josephson's junc- 
tions Hamiltonian. This result is natural because the 
wells and barriers given in Eq. ^ is a Josephson junc- 
tion lattice. 



tb 


42t/4 




— Ji,i+l 


Ji,i+2 




R<p R\ 


40.0 0.070354 0.13905 0.070357 


8.71 x 10" 


'S 


25.1 549 


30.0 


0.16205 


0.32199 


0.16195 


4.70 x 10~ 


-4 


15.5 229 


20.0 


0.41458 


0.81634 


0.41174 


3.10 x 10" 


-3 


8.91 84.5 


10.0 


1.2046 


2.2608 


1.15497 


2.50 x 10" 


-2 


4.72 26.4 


0.0 


3.7826 


5.8779 


3.1027 


1.87 x ID- 


-1 


2.51 7.95 



Table 1: Couplings K {D ^> by DW method and J tJ by 
Gauss approximation at different tb and t w = —30.0. 
We introduce two ratios. is the ratio between the 

maximum of </4°^ at the center of well and its minimum 
at the center of barrier. Another ratio is defined by 
Rx = A 8 /Ay. At higher t b , the saddle point solution 
in the barriers are very small, R^ is very large, the phase 



variation concentrate more in the barriers, so the assump- 
tion is good that the phase variation inside the well is 
ignored. For lower tb, the saddle point solution in the 
barriers is no longer small and the phase variation does 
not favor concentrating in the barriers. The assumption 
ignoring the phase variation inside the well is no longer 
good. Consider the extreme case % = t w , no block can 
be well defined. At higher %, the ratio R^ is very large, 
taking only the first 7 eigenmodes and neglecting other 
modes is a good approximation. For lower tb, the ratio 
R^ becomes small, the approximation to neglect other 
modes becomes bad. 

This approximation is similar to the phase-only ap- 
proximation to simplify the Ginzburg-Landau model to 
XY model 13], in which the modulus of each blocks are 
fixed and only their phases are allowed to fluctuate. 

We summarize the Gauss approximation method as 
follows: 

(1) . Obtaining the saddle point solutions and spatial 
ranges of blocks by DW method. 

(2) . Solving the eigenmodes of <f> y . 

(3) . Defining the block functions as in Eq. (|23p and 
expanding the Hamiltonian as in Eq. (|28|). then we can 
get the couplings. 

C. Gauss approximation: expansion around the 
excited state 

We can also expand the GLW Hamiltonian near the ex- 
cited states with the above method. Consider the excited 
state shown in Fig. (2a). The excited state solution is 
obtained by assign the initial value be negative in the 4th 
barrier and positive at other sites. Therefore we denote 

it by 4> x ^- Let 

4> X = <f> X - , 4 = CPy - ^) = 0y , (30) 

and we consider the Gauss approximation that 

H ^F 4 + 5H X + 5H y (31) 

where F4 is the free energy of excited state (j) x 4 ^ solution 
and 

SH X = d^-[(^-r + it + 3(4 4) ) 2 ]^] (32) 

6H v = £*\[(^f + [t+(<f>W) 2 }4>l], (33) 

where the quartic terms are omitted. 

The eigenmodes of 4>y satisfies the following equation 

- < ^ + [t+{€ ) ?]Vn=KVn (34) 

where A„ is eigenvalues and tp n are the eigenfunctions. 
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FIG. 2: (Color on line) (a) The excited state solution 
(b) The first 7 eigenvalues of Eq. (|34l) . (c) The eigenfunctions 
(pi,(f2- (d) The eigenfunctions <pi,ip5,(p6,<fiT- 



The first 7 eigenvalues are shown in Fig. 2(b) and some 
eigenfunctions are shown in Fig. 2(c) and 2(d). The first 
two eigenvalue are negative. This means that the first 
two modes can cause the free energy to decrease. 

The third eigenfunction ip 3 has eigenvalue A3 = and 
it satisfies that 



14 



(35) 



Table 2: Couplings K {D) by DW method and J l3 by 
Gauss approximation at different t\, and t w — —30.0 for 
the excited state 6 

For B L <C 1 , the effective Hamiltonian of the fluctuation 
about this excited state can be expanded into following 
XY model approximately 



SH„ 



m where 



i<j 



ld G) [(cos(0 Oi + 9i) - (e oj + 9j) - l]/2 (38) 



'Oi — 




(39) 



Comparing Eq. ([28) and ([38]) . we get 



J 3 , 4 = J 4 , 3 = J 4 , 5 = J 5 ,4 - KjQ/4, 

Ji,i+i = Ji + l,i = -K\% x /^ for i^3,4 (40) 

In addition to the difference between Kff\ 1 and K^K, 
the lattice translational invariance is also broken in this 
expansion. However for high tb, the difference between 
two methods becomes very small and the breaking of 
lattice translational invariance also becomes very small. 

We also studied the expansion near other excited 
states, the conclusion is similar. Therefore we show that 
the effective Hamiltonian is given by Eq. (fTT]) with con- 
tinuous order parameter a defined by Eq. (fT2"[) . 



because 4>^ ' also satisfies the saddle point equation (3) 
and Eq. ([34)) with A = 0. This mode corresponds to a 
global rotation. 

Similarly we introduce the block functions 



m) = 



ii 4) ; 2(i - 1) + 0.5 < C < 2i 
0; other cases. 



0.5, 



and assume that 



(36) 



(37) 



1=1 



Here it should be noted that the phase of the 4th block 
is 7r + #4 rather than 8 '4. 



-Jr, 



-Jl: 



J: 



2.:s 



J-3 



40.0 0.070354 0.070357 0.70357 0.070f 76 0.070510 
30.0 0.16205 0.16195 0.16195 0.16097 0.16306 
20.0 0.41458 0.41174 0.41174 0.40510 0.42221 
10.0 1.2046 1.1549 1.1548 1.0965 1.2866 
0.0 3.7826 3.1012 3.0891 2.4233 4.7582 



D. Gauss approximation for a real random 
temperature 

The periodicity in the above discussion is not essential. 
We apply this method to real random temperature cases. 
The couplings obtained by two methods agree with well 
for weak couplings. 

We show a typical example in Fig. 2. The ground 
state saddle point solution is shown in Fig. (2a). There 
are 18 blocks. After solving the Eq. p7[) with the saddle 
point solution shown in Fig. (2a), we get the eigenvalues 
and the eigenfunctions. We divided the systems into 18 
blocks and denned functions similar to Eq. (1231) . Then 



we expand the Hamiltonian of the <j> y fluctuation and 

(G) 

get the couplings K\a ■ The comparison between the 18 
couplings obtained by the two methods are given Fig. 
(2b). 

The couplings obtained by two methods agree well in 
the range from 10~ 4 to 1. The agreement is not good 
for K 5 £, Kq i, K'j g and K w ll because the saddle point 
solution is not small in the regions between these couples 
of blocks. 
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2 4 6 8 10 12 14 16 18 20 
i 



FIG. 3: (Color on line) (a) The ground state solution for 
certain random temperature realization, (b) The couplings 
between the adjoined blocks obtained by two methods. The 
black squares represent K^f^-y and the red circles represent 



IV. TWO DIMENSIONAL CASES 

We consider the well lattice which is defined by 

{ 21 - 1 < £ < 21 
" 2m-l<7]<2m (41) 
ti,; other cases. 

where I, m = 1, 2, 3 and i = (m — 1) * 3 + 1. The subscript 
off, is the label of the well, i.e., for convenience, we label 
the 9 wells with i — 1,2, ■ ■ ■ , 9 as shown in Fig. 4(a). 

Assume the solution along x-direction, the saddle point 
equation is given by 

-^ + ^+^,^ + ^=0. (42) 

where the usual periodic condition is used. 

Expanding the Hamiltonian about this saddle point 
solution, the eigenmodes of <p y satisfy 

-^ + ^ + [*+(^ 0) ) 2 K=An^ (43) 

where A„ is eigenvalues and ip n are the eigenfunctions. 
Through the similar method given in section II, we obtain 
the couplings in Gauss approximation. 

We first consider the simple case with uniform t\ — 
t w = —30.0 for i = 1, 2, 3, • • • , 9. The ground state so- 
lution for this well lattice is shown in Fig. (4b). It can 




FIG. 4: (Color on line) (a) Well lattice, (b) The ground state 
solution, (c) The eigenfunction tp^. (d) The eigenfunction 
ifie- The temperatures are tt = 30.0, t w = —30.0 for (b,c,d). 



be seen that there are 9 blocks, which are the 9 cells 
of the well lattice. Using DW method, we can get the 
couplings between these blocks. Due to the lattice trans- 
lational invariance, we can use Ky 2 and Ky 5 represent 
the coupling between nearest neighbors and next nearest 
neighbors respectively. 



t b Ktf/A -Ji, 2 Kl%>/2 -Ji, B R x 
40.0 0.035649 0.035720 0.000234 0.000199 187 491 
30.0 0.086821 0.087172 0.00104 0.000847 96.5 203 
20.0 0.24358 0.24549 0.00591 0.00455 34.6 72.3 
10.0 0.84624 0.85015 0.0461 0.0336 11.2 19.9 
0.0 3.6426 3.2314 0.348 0.253 3.28 4.75 



Table 3: Couplings K (D *> by DW method and Jy by 
Gauss approximation at different t\, and t w — —30.0. 

As shown in Fig. (4c), for the eigenfunction ip$, the 
phase variations in the wells are much smaller than in 
the barriers. The first 9 eigenfunctions have this feature 
and other eigenfunctions do not have. As a example, 
we show ipie in Fig. (4d). And for tb > 0, the first 
9 eigenvalues are much smaller than other eigenvalues. 
Therefore we take the first 9 eigenmodes into account 
and obtain the couplings between the clocks as we do 
in the preceding sections. The comparison between DW 
method and Gauss approximation are given in Table 3 
for different £(,. The data in Fig. (3b, 3c, 3d) are for 
t b = 30.0. 
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tb 


10.0 


20.0 


30.0 


nian with random temperature. The final conclusion 


ij 


<>/4 - 


m ] k^ia - 


47 K%>,4 -J%> 


that 



12 0.6062 

23 0.8359 

31 1.1650 

45 0.6208 

56 0.7767 

64 1.5560 

78 0.8220 

89 0.3816 

97 0.4866 

14 1.0230 

47 0.3634 

71 1.7770 

25 1.1410 

58 0.4452 

82 0.6657 

36 0.9424 

69 0.5280 

93 0.6061 



0.6384 
0.8466 
1.1270 
0.6649 
0.8107 
1.4610 
0.8111 
0.4314 
0.5384 
1.0130 
0.4064 
1.6400 
1.1080 
0.4790 
0.7261 
0.9379 
0.5488 
0.6679 



0.1610 
0.2400 
0.3630 
0.1651 
0.2189 
0.5194 
0.2343 
0.0899 
0.1217 
0.3087 
0.0861 
0.6113 
0.3522 
0.1091 
0.1794 
0.2778 
0.1358 
0.1591 



0.1651 
0.2429 
0.3617 
0.1711 
0.2246 
0.5114 
0.2345 
0.0955 
0.1286 
0.3099 
0.0903 
0.5977 
0.3509 
0.1136 
0.1886 
0.2792 
0.1384 
0.1679 



0.0545 
0.0854 
0.1358 
0.0557 
0.0766 
0.2031 
0.0831 
0.0278 
0.0390 
0.1132 
0.0269 
0.2440 
0.1312 
0.0353 
0.0604 
0.1005 
0.0450 
0.0527 



0.0552 
0.0859 
0.1359 
0.0566 
0.0776 
0.2022 
0.0832 
0.0286 
0.0401 
0.1135 
0.0264 
0.2421 
0.1311 
0.0359 
0.0619 
0.1009 
0.0454 
0.0541 



Table 4: Couplings K\f > by DW method and J^> by 
Gauss approximation with nonuniform ij at different th- 
in fact only the inhomogeneity rather than the peri- 
odicity of the temperature is essential. We also con- 
sider the random temperature cases with ti,t%, ■ ■ ■ ,tg be- 
ing -30, -20, -50, -40, -17, -60, -35, -25, -15 respec- 
tively. There are 18 couples of nearest neighbored blocks. 
In the table 4, we present the couplings between nearest 
neighbored couples at % = 30.0, 20.0, 10.0. As we can see 
that for tb — 30.0 the differences between the couplings 
obtained by two methods are much smaller than those at 
t b = 10.0. 



(G) 



V. SUMMARY 

The continuous fluctuation about the saddle point so- 
lution is studied for XY-type Ginzburg-Landau Hamilto- 



(1) The effective Hamiltonian of blocks is a XY model. 

(2) The DW method provides the excited states with 
block's phase being only 0, n, and the Gauss approxima- 
tion can describe the continuous phase fluctuation near 
the saddle point solutions. 

(3) The couplings obtained by these two methods agree 
with each other well for the weak coupling cases. 

DW method has a great advantage comparing with the 
Gauss approximation. For DW method, the size of the 
grid in the numerical calculation can be as large as 2000 x 
2000, while in the Gauss approximation, the size of grid 
can only be as large as 120 x 120 and the computing time 
is very long because of diagonalizing matrix. Therefore 
using DW method, the couplings between blocks can be 
conveniently calculated and the statistical properties of 
the couplings can be studied. 

Our conclusion is consistent with the recent experi- 
ments and theoretical studies. Appearance of granular 
structures self-organized in homogeneously-disordered 
SC is discovered by the experiment [bH ] and shown 
by theoretical studies near the quantum critical point 

Recently the excited state solutions of Bogliubov- 
de Gennes equations are solved for two dimensional 
neg ative-U Hubbard Hamiltonian with on-site disorder 
}17| . The excited states show that the system is self- 
organized into blocks. DW method is used to obtain 
the couplings between blocks. The authors claimed that 
the effective Hamiltonian between these blocks should be 
XY-type. The argument in this paper can be regarded 
as its corroborative evidence. 
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